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The lifetime measurement of n + n~ atoms (Ai%) with 10% precision provides in a model independent way the difference between the 
S-wave Tin scattering lengths for isospin and 2, \uq — a 2 \, with 5% accuracy. The scattering lengths gq and a 2 have been calculated in 
Chiral Perturbation Theory (ChPT) with a precision better than 2.5%. Therefore, such a measurement will be a sensitive check of the 
understanding of chiral symmetry breaking in QCD by giving an indication about the value of the quark condensate, an order parameter 
of QCD. The method of observation and lifetime measurement of A 2% is discussed. Data on A 2x production in Pt, Ni and Ti targets and 
preliminary results on the lifetime measurement are presented. 
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Introduction 



Pionium or A 2n is a hydrogen-like atom consisting of n + 
and n~ mesons. The lifetime of this atom is inversely pro- 
portional to the squared difference between the S-wave %% 
scattering lengths for isospin and 2, |oq — a%\. This value 
is predicted by chiral perturbation theory (ChPT), and a 
measurement of the % + n~ atom lifetime provides a possi- 
bility to check predictions of ChPT in a model-independent 
way. 



1 Present status of nn scattering length in- 
vestigation. 

In chiral perturbation theory (ChPT) the 2-loop calculation, 
exploiting Roy equations, leads to a precise prediction for 
\a Q -a 2 \ at the level of 1.5%: \a Q -a 2 \ = 0.265 ±0.004 [ 



Experimental data on tztz scattering lengths can be obtained 
in an indirect way via the following processes: 



1. K e 4 decay, K + =>• 7C + 7C~e + V e , using a model- 
independent analysis. The most resent result ob- 
tained from 4 x 10 5 K e 4 decays have been col- 
lected in the experiment E865 (BNL). The analysis 
yields the S-wave nn scattering length «o = 0.216 ± 
0.013(stat) ±0.002(syst) ±0.002(theor) [0. 

2. nN — * nnN reaction. The latest measurement gives 
a value a = 0.204 ± 0.014(stat) ± 0.008(syst) [0. 
In this case the experimental procedure is not strictly 
unambiguous, since the contribution of non-OPE 
components (A, TV*) can effect the reliability of the 
result [j3J. 



2 Theoretical motivation 



The n + n~ atom decays by strong interaction mainly into 
7r°7r°. The branching ratio of the alternative decay mode 
A-2n —> 2y is at the level of 4 • 10~ 3 . There is a relation 
[IS Hi between r 27r o and nn scattering lengths: F 2n o = 
C ■ \ao — a2\ 2 ■ 

ChPT at next-to-leading order in isospin breaking pro- 
vides correction to decay width I™q = r 2jr o(l + 5r), 5r = 
(5.8 ± 1.2)% . Using this correction the lifetime of n + 7t^ 
atoms is predicted to be [|6): 



X= (2.9±0.1)-10~ 13 s 

A measurement of the atom lifetime with an accuracy of 
10% allows to determine the corresponding n% scattering 
length difference with 5% accuracy. 

3 Method of lifetime measurement 

The A 2 n are produced by Coulomb interaction in the fi- 
nal state of n + n~ pairs generated in proton-target interac- 
tions [ElEJ. After production A 2n travel through the target 
and some of them are broken up due to their interaction 
with matter: "atomic pairs" are produced, characterized by 
small pair cm. relative momenta Q < 3 MeV/c. These 
pairs are detected in the DIRAC setup. Other atoms annihi- 
late into n° n°. The amount of broken up atoms iia depends 
on the lifetime which defines the decay rate. Therefore, the 
breakup probability is a function of the A 2k lifetime. 

The dependence of F\, r on the lifetime T is determined by 
the solution of differential transport equations [0. In 
Fig-Hthe lifetime dependence of Pj,r is presented for three 
different targets used in the DIRAC experiment. The nickel 
target provides the best statistical accuracy for the same 
running time. 
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Figure 1. Dependence of the breakup probability f\, r on Ain life- 
time for three targets used in the DIRAC experiment: platinum of 
26 fim, nickel of 94 flm and titanium of 247 flm thickness. 
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Figure 2. Simulated distributions of "atomic pairs" (solid line), 
"Coulomb pairs" (dashed line) and "non-Coulomb pairs" (dotted 
line). 



Also pairs from short-lived sources are generated in 

free state. Such pairs ("Coulomb pairs") are affected by 
Coulomb interaction, too. The number of produced atoms 
(Na) is proportional to the number of "Coulomb pairs" 
(N c ) with low relative momenta (Na — K -N c ). The coef- 
ficient K is precisely calculable. And there are n + n~ pairs 
from long-lived sources (rj,K^,...). Such pairs, not af- 
fected by final state interaction, are named "non-Coulomb 
pairs". 

In Fig. 12 simulated distributions of "atomic pairs", 
"Coulomb" and "non-Coulomb" pairs are presented. For 
relative momentum Q > 3 MeV/c "atomic pairs" are prac- 
tically absent, but in the region Q < 2 MeV/c a fraction of 
them is essential, and this region provides the best effect- 
to-error ratio. 

The aim of DIRAC is to measure the Ai % breakup proba- 
bility Pbr(f). f\}i(t) is the ratio between the observed num- 
ber of "atomic pairs" and the number of produced n + n~ 
atoms which is calculated from the measured number of 
"Coulomb pairs". 

4 Experimental setup 

The purpose of the DIRAC setup (Fig.|3J is to detect n + n~ 
pairs with small relative momenta F II 01 . This setup is lo- 
cated at the CERN T8 beam area (East Hall). It became 
operational at the end of 1998 and uses the 24 GeV proton 
beam from PS accelerator. 

In Tabled the main setup features are summarized. The 
setup resolution is shown for the projections of the relative 
cm. momentum Q: Qi is a projection along the total pair 
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Figure 3. DIRAC setup. MSGC are microstrip gas chambers, 
SFD is a scintillating fiber detector and IH is a scintillation ion- 
ization hodoscope. Downstream the spectrometer magnet there 
are two identical arms Tl and T2. Each arm consists of drift 
chambers (DC), vertical (VH) and horizontal (HH) scintillation 
hodoscopes, threshold Cherenkov counters (CH), shower detec- 
tors (PSH) and scintillation muon detectors (MU) 

momentum in the laboratory system, Qx and Qy are the X- 
and Y- projections of the transverse component of Q. 



5 Experimental data 

The setup detects pairs of oppositely charged particles. 
Fig-E]shows the event distribution over time difference be- 
tween positive and negative particles. The peak in the cen- 
tral part contains pairs of particles which are generated in 
the same proton-nucleus interaction ("real pairs"). These 
pairs consist of "atomic", "Coulomb" and "non-Coulomb" 
pairs. The wings of the distribution are formed by "acci- 
dental pairs" consisting of particles from different proton- 
nucleus interactions and so not affected by final state in- 
teraction. Therefore, the g-distribution of "non-Coulomb" 
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Table 1. Setup features 



Angle of channel to proton beam 


5.7° 


Channel aperture 


1.2 KT 3 sr 


Magnet bending power 


2.3 T m 


Momentum range 


1-2 <p n <7GeV/c 


Resolution on pair relative 


Oq x = Gq y = 0.4 MeV/c 


momentum projections 


Oq l = 0.6 MeV/c 



and "accidental" pairs are very similar. Only a laboratory 
momentum distribution for "accidental pairs" is more hard 
than one for "real pairs" due to energy conservation. It 
causes the difference for big values of pair relative mo- 
mentum Q. But in DIRAC only events with small relative 
momenta are used. Admixture of "accidental pairs" under 
the central peak is subtracted by using the widths of the 
intervals Afj , At 2 and At 3 . 
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Figure 4. Distribution of events over time difference measured by 
vertical scintillation hodoscope (time of negative particle is sub- 
tracted from time of positive one). Regions Afi and Af3 contain 
pure "accidental pairs" and A/2 one contains a mixture of "real" 
and "accidental" pairs. 



Eq. ^ shows the experimental correlation function R(Q). 
This function consists of two parts. One part is a con- 
stant N ■ f as the result of dividing the distribution of "non- 
Coulomb pairs" (dN^ /dQ) by the distribution of "acci- 
dental pairs". Here N is a normalization parameter and / 
the fraction of "non-Coulomb pairs". 



ilN rF 



HQ) 



dQ 
dQ 



dN NC . 

real 

dQ 



dN c . 

real 

dQ 



dQ 



= N-[f+A c (Q)} 



(1) 



The second part of R(Q) is Q-dependent: the number of 
"Coulomb pairs" (dN^/ dQ) with very small Q is in- 
creased by Coulomb interaction. It should be emphasized 
that the Coulomb factor Ac(Q) (Eq. ^ is smeared com- 
pared to the theoretical one due to multiple scattering in 
the target and finite setup resolution. 



The experimental relative momentum distributions for 
"real" and "accidental" pairs are shown in Fig. [5] These 
data were collected with nickel target in 2001. 




Figure 5. Experimental distributions of K + 7l pairs over relative 
momentum Q for "real" pairs (a) and "accidental pairs" (b). 

Fig-Elshows the experimental correlation function (Eq.[Q. 
In the region Q > 3.5 MeV/c this function was fitted by 
an approximated correlation function which takes into ac- 
count "non-Coulomb pairs" (constant term) as well as 
"Coulomb pairs" (Q-dependent term), the multiple scat- 
tering in the target and also the setup resolution. The 
fraction of "non-Coulomb pairs" / is a free parameter in 
the fit procedure (Eq. There is a good agreement be- 
tween the experimental and approximated correlation func- 
tions for Q > 2.5 MeV/c. On the other hand in the region 
Q < 2.5 MeV/c the experimental points are higher than the 
approximation function describing "free pairs" only. This 
excess is the signal of "atomic pairs". 

To calculate the number of "atomic pairs" the approxi- 
mated correlation function has been transformed to a Q- 
distribution by multiplying it with the distribution of "ac- 
cidental pairs" (Eq. 0. Differences of experimental and 
approximated distributions are shown in Fig.0for data col- 
lected in 1999, 2000, 2001 and in 15 days of 2002. 

Simulation shows that analysis of Qi-distributions pro- 
vides a decreasing of the systematic error by a factor ~ 2 
and a small increasing (30%) of the statistical error in com- 
parison with analysis of ^-distributions. 

The distribution over the absolute value of Qi for a mix- 
ture of "real" and "accidental" pairs and the distribution of 
"accidental pairs" only are shown in Fig. [8] for data col- 
lected with the nickel target in 2001. Data were selected 
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Figure 6. (a) The experimental correlation function of k + k~ 
pairs R(Q) (points with error bars) is fitted by an approximated 
correlation function (smooth curve) in the region Q > 3.5 MeV/c. 
The excess for Q < 3.5 MeV/c occurs due to "atomic pairs", (b) 
The approximated correlation function is the sum of correlation 
functions for "Coulomb" (points) and "non-Coulomb" (dashed 
line) pairs. 

with the criterion on transverse component of relative mo- 
mentum Qt < 3 MeV/c. Therefore a small value of Ql 
means a small value of Q. It is seen that the distribution of 
real pairs shows a peak at Ql = 0. 

The distribution of "real pairs" (Fig|9t) is fitted by the sum 
of approximated distributions of "Coulomb" and "non- 
Coulomb" pairs in the region \Ql\ > 2 MeV/c where 
"atomic pairs" are absent. Approximated distributions are 
made on the base of "accidental pairs" applying a corre- 
lation function which takes into account the Coulomb in- 
teraction in the final state, multiple scattering in the tar- 
get, the resolution of the setup and difference of labora- 
tory momentum spectra of "real" and "accidental" pairs. 
The fraction of "non-Coulomb pairs" and the normaliza- 
tion factor are free fit parameters. There is a good agree- 
ment between the experimental and approximated distri- 
butions for Ql > 2 MeV/c. On the other hand in the re- 
gion Ql < 2 Me V/c the experimental points are higher than 
the values of the approximation function describing "free 
pairs" only. This excess provides a number of "atomic 
pairs" and is shown in Fig |5p separately. The number of 
"Coulomb pairs" is calculated in the same fit procedure. 

The number of "atomic pairs" ha — 5224 ± 397 together 
with a number of "Coulomb pairs" and the theoretically 
known coefficient K allows to estimate a breakup proba- 
bility value of .Pbr( T ) = 0.462 with a statistical error 9.2%. 



Figure 7. Signals of "atomic pairs" in the region Q < 3 Me V/c 
are obtained with different targets. For 2002 only data collected 
in 15 days are analyzed. 

There is also a systematic error due to accuracy in descrip- 
tion of multiple scattering in the target, membranes and 
detector planes (precision needed is 1 %) as well as in de- 
scription of the detectors response. 

Taking into account the dependence of breakup probability 
on lifetime (Fig.llOt a preliminary estimation of the 7t + n~ 
atom lifetime yields 

T = [3A2t™l(stat) ±l.(syst)} ■ 10~ 15 s 
Conclusions 

The DIRAC collaboration has built up a double arm spec- 
trometer providing a resolution of 1 Me V/c for low relative 
(Q < 30 Me V/c) momentum of particle pairs and has suc- 
cessfully demonstrated its capability to detect n + n~ after 
2 years of running time. 

In order to decrease systematic errors, dedicated measure- 
ments with a multi-layer nickel target and measurements of 
multiple scattering in all detectors and setup elements have 
been performed at the end 2002 and during present run of 
2003. 

Preliminary results have been achieved by analyzing data 
collected in 2001. The statistical accuracy in the lifetime 
determination reaches 25% and the systematic one is 30%. 
The analysis of data collected in 2000 and 2002 together 
with systematic error reduction allows us to improve the 
accuracy up to the level of 14%. 
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Figure 8. Experimental distributions of k + k~ pairs over Q L - 
projection of pair relative momentum for mixture of "real" and 
"accidental" pairs (a) and "accidental pairs" (b). 
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Figure 10. Preliminary estimation of k + k atom lifetime with 
measured breakup probability. 
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Figure 9. Experimental and approximated (dashed line) distri- 
butions of "real" k + k~ pairs over Qi (a) and signal of "atomic 
pairs" in the region <2l < 2 MeV/c is obtained using data collected 
in 2001 with a nickel target (b). 
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